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The transmission (T) and reflection {R) coefficients are studied in periodic systems 
and random systems with gain. For both the periodic electronic tight-binding model 
and the periodic classical many-layered model, we obtain numerically and theoreti- 
cally the dependence of T and R. The critical length of periodic system L^, above 
which T decreases with size of the system L while R approaches a constant value, is 
obtained to be inversely proportional to the imaginary part e" of the dielectric func- 
tion e. For the random system, T and R also show a nonmonotonic behavior versus 
L. For short systems (L < Lc) with gain < InT >= {Ig^ — ^q^)L. For large systems 
(L >> Lc) with gain < InT >= —{^g^ + ^q"^)!^- Lci Ig and £^re the critical, gain 
and localization lengths, respectively. The dependence of the critical length Lc on 
e" and disorder strength W are also given. Finally, the probability distribution of 
the reflection R for random systems with gain is also examined. Some new very 
interesting behaviors are observed. 
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I. INTRODUCTION 



While the study of locahzation of classical and quantum waves in random disordered 
media has been well understood flj-Q , recently, the wave propagation in amplifying random 
media has been pursued intensively p|-[T3[l. Some interesting results have been predicted, 
such as, the localization length of a random medium with gain |^, the sharpness of back 
scattering coherent peak |^,rU|, the dual symmetry of absorption and amplification [^], the 
critical size of the system and the probability distribution of reflection ||^. Numerically, 
two kinds of models are studied, one is the electronic tight binding model [^,^, the other 
is the many-layered model of classical waves [§] . Theoretically, a lot of methods are used 
to get these results, such as the diffusion theory and the transmission matrix method 

0. Most of these studies are for homogeneously random systems which are generated by 
introducing the disorder into the continuous system, and the medium parameter, such as 
the dielectric constant, is assumed to vary in a continuous way ||5|J^. But the periodically 
correlated random systems which are generated by introducing the disorder into a periodic 
system, such as a photonic-band-structure , have not been studied adequately. 

With gain, will such random systems with periodic background behave similar as the 
homogeneously random system? Both experimentally and theoretically, the study of such 
system is very important in understanding the propagation of light in random media. These 
type of photonic-band-structure systems are widely used in experiments Theoretically, 
just as S.John P| argued, the localization of a photon is from a subtle interplay between 
order and disorder. For the periodically correlated random systems with gain, the periodic 
background plays the order role, and now its interplay with not only disorder but also with 
gain should be a very interesting new topic. 

In this paper we address both the electronic tight binding model and the many-layered 
model of classical waves. We first compare the numerical results of periodic amplifying 
system with what we can predict theoretically by the transfer matrix method. It's surprising 
to get most of the universal properties, such as critical length and exponential decay of 
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transmission, of homogeneously random system ||^ from a periodic system too. 

With the help of some theoretical arguments and numerical results, we suggest that the 
length ^1 = \1/ Im{K)\, where K is the Bloch vector in a periodic system with gain, to 
replace the gain length Ig = |l//m(fc)| introduced in p[. This is more reasonable since the 
correlated scatterers in a periodic system can make the paths of wave propagation much 
longer in the system. We also think that it is actually the Bloch wave instead of the 
plane wave which propagates in the system. Then we introduce disorder into these periodic 
systems and calculate their properties. Our numerical simulations for both models show that 
periodically correlated random systems give similar behaviors as that of the homogeneously 
random systems studied previously. But in some cases, we get new interesting results for 
the localization length ^, the critical length Lc and the probability distribution P{R) of 
reflection. All these new results are related to the periodic background of such systems. 
We also examine the results of the transmission coefficient T for short {L < Lc) systems. 
Our numerical results show that the formula of the transmission coefficient of media with 
absorption can be generalized to the transmission coefficient of short systems with gain, if 
we replace the gain length Ig (or ^i) with the negative of the absorption length —la in the 
formula. To explain our new results of the critical length L^, we compare the two basic 
theories for obtaining the critical length, the Letokhov theory [^] and the Lamb theory ||16 



and we get some new theoretical results of critical length which are in good agreement with 
our numerical results. The behavior of the distribution of the probability of reflection P{R) 
is much more complex than the theoretical prediction of homogeneously random system . 
We flnd that the periodic background influences strongly the general behavior of P{R)- 

The paper is organized as follows. In Sec. II we introduce the two theoretical models we 
are studying. The results for the periodic systems with gain are presented in Sec. Ill, while 
in Sec. IV the results for the random systems with gain are given. Also in Sec. IV we present 
our new theoretical and numerical results for the critical length Lc. In Sec.V the results for 
the probability distribution of reflection coefficient R for both models are presented. Finally, 
Sec. VI is devoted to a discussion of our results and give some conclusions. 



II. THEORETICAL MODELS 



A. Many-layered model of classical wave 

Our periodic many-layered model of classical wave consists of two types of layers with 
dielectric constant ei = Eq — ie" and €2 = X^o "~ i^" and thicknesses a = 95nm and 60 = 
120nm respectively, where the negative part of dielectric constant, i.e. e" > 0, denotes the 
homogeneous amplification of the field. We have tried a lot values for x, such as 1.5, 2, 3, 5, 
6, and get no essential difference in our results for different values. In this paper we choose 
X = 2, i.e. Re{ei) = Eq = 1 and Re{e2) = 250 = 2. The system has L cells. Each cell is 
composed by two layers with dielectric constant £1 and 62 respectively. Without gain, we 
obtain that the wavelength range of the second band of this periodic system is from 247 
nm to 482.6 nm ( the first band has a range from 592 nm to infinite). So we choose the 
wavelength 360 nm to represent band center, the wavelength 420nm as a general case, and 
the wavelength 470 nm to represent the band edge. 

To introduce disorder, we choose the width of second layer of the nth cell to be random 
variable 6„ = 6o(l + W^7); where W describes the strength of randomness and 7 is a random 
number between (—0.5,0.5). The whole system is embedded in a homogeneous infinite 
material with dielectric constant equal to ^o- 

For the ID case, the time-independent Maxwell equation can be written as: 

Suppose that in the medium with dielectric constant Si and the medium with dielectric 
constant £2, the electric field [H is given by the following expressions: 

E2n{z) = C^e'''^'-'-^ + Dr.e-'''^'-'-^ (2) 

Using the appropriate boundary condition (continuity of the electric field E and of the 
derivative of E at the interface), we obtain that: 
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where k — ^ ~ ^V^- 

Prom the product of these matrices, M(L) = Ilf -^n, we can obtain the transmission and 
reflection ampUtudes of the sample, t{L) = and r{L) = For each set of parameters 
{L,W,e"), the reflection coefficient R — |rp and the transmission coefficient T — \t\'^ are 
obtained from a large number of random configurations. We have used 10,000 configurations 
to calculate the different average values of R and T, and 1,000,000 configurations to obtain 
P{R)- Our numerical results show that the localization length for a system without gain 
behaves oc for this model, and are in agreement with previous workers. 



B. Electronic tight-binding model 



For the electronic tight-binding model, the wave equation can be written as: 
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where 
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(6) 



e„ = Wj—ir), where W describes the strength of randomness, 7 is a random number between 
(—0.5, 0.5) , r] > corresponds to amplification and 0„ is the wave function at site n. The 
length L of the system is the total lattice number of the system. The system is embedded 
in two identical semi-infinite perfect leads on either side. For the left and the right sides, we 
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have (f)o = l + r{L) and (pL+i = t{L)e'^^^+^\ We can obtain reflection amplitude r(L) and 
transmission amplitude t{L) by the products of matrices, M{L) = H^Mn. 
t(L) = -2tstn{k) 

""^ ^ Mne-^'= + Mi2-M2i-M22e^^'' ^' 

where k = arccos{^). 

When W = and without gain, the model is a periodic one with only one band spanning 
in energy between -2 and 2. Notice that the hopping matrix elements in Eq.(^ are equal to 
one, which is our unit of energy. So we choose = to represent band center, = 1 as a 
general case, = 1.8 to represent band edge. 

Similar as the many-layered model, for each set of parameters (L, W,ri), 10,000 random 
configurations were used to obtain a average value of R and T, and one million random 
configurations for P{R)- Theoretical and numerical results give that the localization length 
for a system without gain behaves ^0 oc 1/W'^, in agreement with previous workers. 

III. PERIODIC SYSTEMS 

Almost all the properties of the periodic systems of both the many-layered model and 
the tight-binding model can be predicted theoretically. 

A. Classical many-layered model 

For long systems (L ^ of the many-layered model we have that: 

L->oo oL 

= 2Im{K) oc e" (8) 

where K is the Bloch vector which is a complex number now, and satisfies cosK = 
cos{ka)cos{qbo) — -|- ^)sin{ka)sin{qbo). Because Im{K) < 0, the transmission coeffi- 
cient T is decaying exponentially for a long system. 
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For a short system, L < L|?, we have: 

d\nT 



l/C,c^2\C\\Im{K)\ 



dL 

^2\Im{K)\ (9) 

Where C = —[sin{ka)cos{qbQ) + cos{ka)sin{qbo)]/[2sin{K)], and \C\ is larger than but very 
close to 1 when wavelength is at the band center, and become bigger when wavelength 
approached the band edge. 

So the slope of InT vs L for a short periodic system is almost same as the negative 
value of the slope for the long system. The slopes of In T at both sides of the maximum are 
approximately symmetric. In Fig. la, we can see that, when L < T increases vs L with 
the slope and get to a maximum at L° and decays exponentially when L > with 
the slope 

From the behavior of the theoretical expressions of T or i? , when T or R goes to infinite, 
we can obtain analytically that is given as: 

i2==6.n(^) (10) 

where C is same as defined above in Eq.(^, and \C\ is close to one. From the property of 
\C\ discussed above, we can see that L° >> at the band center , and becomes smaller 
when wavelength approaches the band edge. The value of \C\ is almost independent of gain, 
so is parallel to l/e" or ^i. We have shown that L^e" is almost a constant for a given 
wavelength, and our numerical results agree very well with the theoretical prediction. 

The reflection coefficient gets to a maximum value at too, and fluctuates a lot with 
the size L of the system. When L approaches infinity, R reaches a saturated value. The 
saturated value of R is given by: 

hm R = RoC^ (11) 

\stn{K){\C\ - l)f ^ ' 

So Ro is almost independent of gain or ^i. Fig. 2 shows that indeed R increase when L < L'^, 
gets to its maximum value and fluctuate violently at L^, then approaches a saturated value 
which is almost independent of gain. 



B. Electronic tight-binding model 



For the electronic tight-binding model, when E — 0, the In T for the long system can be 
obtained by the use of Eq.(7): 

Similarly as the classical many-layered model, for short system of tight-binding model 
we have: 

ainT 



^ r] (13) 

So the slope symmetry of InT at both sides of L° still exists. In Fig. lb, we can see a 
similar behavior as in Fig. la, when L < L°, InT change vs L with the slope of and 
gets maximum at L°, then it begin to decay exponentially with a slope of l/^i. 

Assuming that the theoretical expression of T, given by Eq.(7), is infinite, we can obtain 
that L° is given by: 

L° ~ ^(In4-ln7?) ~ 2^iln(4^i) (14) 

We also shown that L'^rj + 2hirj ys rj, for £" = 0, is a constant for different gain and indeed 
find out that the theoretical prediction given by Eq.(14) agree very well with the numerical 
results. 

The reflection coefficient R approaches a saturated value as L goes to infinite, but the 
saturated value of Rq is not a constant independent of the gain as in the case of classical 
many-layered model. This is clearly seen in Fig. 3 where we plot Ini? vs L. Notice that 
the Ini? curves increase vs L when L < L^, get to a maximum at L°, and then approach a 
saturated value when L goes to infinity. Similar results were obtained for E ^ 0. 



IV. RANDOM SYSTEMS 



In Fig. 4, we give the general behavior of average value < InT > vs L for both models. 
We can see the different behaviors for L < Lc and L > Lc- When l/^i > l/^o and L < Lc, 
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< InT > increase vs L from origin with a slope which is defined as 1/^' and when L > Lc, 

< InT > decrease vs L with a slope —1/^. But when < l/.^0) < > "will decrease 
monotonically, at first with the slope l/^i = — l/|^i|, at Lc, there are a turning point and 
slope changes to —1/^. We will study the values of ^' , ^ and Lc in this section. 

It was first suggested by Zhang |^ that the localization length ^ of a long random system 

with gain will become smaller than the localization length of the random system without 

gain. In particular he suggested that: 

1 dlnT 1 1 ,^ . 

- = hm — — = ^ + ^ 15 

where 6 is the localization length of the system without gain, Ig is replaced by 6 
original formula of Zhang because of the periodic background of our systems. 

We have numerically calculated 1/.^ for different cases of disorder, gain and fre- 
quency (energy), and compare it with 1/6 + shown in Fig. 5a and 5b for the many- 
layers model and the tight-binding model respectively. For most of the cases. Eg. ([T5|) is a 
very good formula. Only when wavelength is on band edge and when both gain and ran- 
domness are very strong, we can see that the numerical results deviate from the theoretical 
prediction, which is the solid line in both Fig. 5a and 5b. 

For short system(L < Lc), the behavior of < InT > vs L is quite different from that of 



long system as shown in Fig. 4. Freilikher et.al. and Rammal and Doucot [17| obtained that 



the transmission coefficient of a random system with absorption is given by: 

<lnT>=(-i-i)T (16) 

where la is the absorption length and 6 is the localization length. For a medium with 
gain, can we just substitute the —la^^ with Ig^^ in the Eq.(|16[) to get following equation? 

<lnT> 1 ,1 1, 

Because of the periodic background of our models, we use 6 to replace Ig in our calculations. 

So far there is no independent verification for this conclusion. After substituting 6 for 
Ig, our numerical results show that Eq.(p!7|) is correct for short systems with gain for both 
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models. When the strength of the disorder is a constant , so is a constant, according to 
Eg. (p!7D , — \/^' should be equal to l/,^o and be constant as the gain varies. We have 
checked this prediction and find indeed that the numerical values are almost the same as 
the ones predicted theoretically. 

From Eq.(|I7p, we can predict the basic features of the length dependence of < InT > 
shown in Fig. 4. When l/^i > l/^o, < InT > will increase with L, and will reach a maximum 
value when L gets to Lc- But if l/^i < l/^o, the < InT > will decrease monotonically, at 
first with a slope of — |l/^i — l/^o| from the origin, at the curve has a turning point and 
the slope changes to — |l/^i + l/,^o|- If — V^o, the curve is almost horizontal for small 
L and begins to decrease with a slope —1/^ at the critical length. This behavior is exactly 
shown in Fig. 4. 

The critical length L^. is one of the most important parameters of a random system with 
gain. For a random system, one of the most important theories is the Letokhov theory 
0. Zhang P generalized the theory and used the no-gain localization length to replace 
the diffusion coefficient D in the Letokhov theory and obtain that the critical length — 



^oig, so that we can clearly see localization effects in the system. But as shown above, 
there is a finite critical length in periodic system when the no-gain localization length 
goes to infinite, so there must be other mechanisms for determining the critical length in 
those systems. We find that when the localization effect is strong enough so that the no- 
gain localization length ^ (-^c)^/'Ci) then the results of the Letokhov theory are quite 
good. But when the system randomness is weak so that is larger than {L^^/^i, then the 
Letokhov theory results are not correct, and we have to use other theories, such as the Lamb 
theory [|16|, which is well known is laser physics, to determine Lc. Next we will compare 



the Letokhov theory with the Lamb theory, and find the expressions of the critical length 
in different cases. 

According to the Letokhov theory P,pn| the field in the system satisfies: 



^^^'^'^ -DVW,t) + "-^ (18) 



dt Ig 
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where D is the diffusion coefficient, c is the speed of the wave. 

Considering the relaxation after long time, the solution [|1^] of Eq.(|18|) is : 

(p{r, t) oc e '9 
^ = -fl(|)V,<) + ^ (19) 

When L = Lc = tt^J ^ , the system is at a critical point. If L < the field will decay 
vs time, but ii L > Lc then the field in the system will become stronger and stronger with 
time. 

We can clearly see that the physical meaning of Lc is the balance point of the gain 
and loss in the system. When the L is less than L^, the photon escaping rate, which is 
determined by is larger than the photon generating rate, which is determined by f- of 
the system, so the photons generated by the stimulated emission can escape from the system 
instantaneously and the system can get to the static state after a long time. If L is larger 



than Lc, gain is larger than loss, and photons will be accumulated in the system [0. Based 



on the Letokhov theory and the weak localization theory |1|,0] results, Zhang p[ generalized 



the Lc to be Lc — y^o^g since D = |/c, where / is the mean free path and ^0 = (2 ~ 4)/. 

In our models, considering the periodic background, we substitute ^1 for Ig first. But 
when the disorder becomes weaker and weaker, the system become almost periodic, ^0 goes 
to infinite, Lc goes to L° instead to infinite. How one can explain this behavior of L^? The 
Lamb theory can give a theoretical explanation of it. In the Lamb theory, a phenomenological 
parameter Q(L), the quality factor which generally is a function of system length L, is 
introduced to show the energy loss rate of the system (also can be thought as the photons 
loss rate of the system). In the Lamb theory, the magnitude of the electric field in a linear 
medium satisfies the following equation: 

At the critical condition, the gain term and loss term are equal. We have: 

= - (21) 
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If we compare this Eq.(pOD with the solution of Letokhov theory Eq.(19), we can find the 
similarity between them. This similarity is from the same physical principle, the interplay 
of loss and gain in the system. From Eq.(^) we can see that the gain term is same as 
the one given by the Letokhov theory, the only difference is from loss term. Generally, ^ 
is a function of the system length, e.g. for Fabry- Perot interferometer ^ oc For 
periodic system Q = Qp, we have oc j too. From the balance of gain and loss, we can 
get that oc in agreement with our results presented in section III. This means that 
in a periodic system the rate of loss is not infinite, although the no-gain localization length 
goes to infinite. The rate is determined by the Qp of the system and we can get a finite L° 
correspondingly. From the L° obtained above and the critical condition given by Eq . , 
we have that the quality factor of the periodic system is given by: 

Q.-^L (22) 

which is independent of no-gain localization length ^q. 

For a random system, things are a little more difficult. The theory of Letokhov doesn't 
give the detailed information of localized modes but it gives a localization related quantity 
D, the diffusion coefficient. According to the localization theory, D is directly related with 
localization length .^o; just as Zhang discussed [§]. Based on the correct results of the 
Letokhov theory in the strong localization case, we can assume that when the disorder is 
strong enough ^ -^c^/^i; the localization effects will dominate the escape rate of photons 
of the system (Our numerical results shown in Fig. 11 support this assumption). Lamb 
theory gives that the Q of a strong random system is determined by the localization effect. 
By comparing the corresponding terms in Eq.(19) and Eq.(EO]), we obtain that: 



where the subscript / is for localized modes, a is a constant of the order of unity and depends 
on the ratio of D and according to the localization theory. For both of the models studied 
here, we find the a can be chosen to be equal to 0.7. From this we can get that the critical 
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length Lc = y-^^o^i — V^o^i which is consistent with the Letokhov theory. Eq.(^3|) is a very 
interesting resuh for laser physics because it is obtained by the comparison of the Letokhov 
and Lamb theories, and it directly gives the relationship of the quality factor Q of a random 
system with the no-gain localization length of the system. 

In the weak disorder limit ^ -^c^/^i; Q ^ Qp and Lc L^. In strong disorder limit 
^0 ^ L^c lii^ Q Qi and Lc V^o^i- For cases where is comparable to both the 

effects of periodic background and randomness will be important to determine the quality 
factor of such a system. Considering the Q as the photon-resistance in the system, and if 
we assume that both effects are independent with each other, we have that the total quality 
factor of the system to be: 

From the critical condition Eq.(^), we have that: 

eo6 



2aL0^'^[2aL0j 

In Fig. 6a and 6b, we compare the theoretical predictions given by Eq.(p5D and by Zhang 
PI with our numerically calculated results for the classical many-layered model and the 
electronic tight-binding model respectively. Our numerical results shown in Fig. 6a and 6b 
strongly support Eq.(P5D to be the correct expression of the critical length Lc for both the 
weak and the strong random limits. In some other cases, the deviation can be as large as 
fifteen percent which is still very good considering the many approximations that have been 
introduced in derivation of Eg . (p5|) . One explanation for this deviation is that the two effects 
that were added in Eq.(|^) are not totally independent, because the correlated scattering 
of the periodic background will affect both Qp and Qi. 

V. PROBABILITY DISTRIBUTION OF REFLECTION COEFFICIENT 



(Ml] 



Pradhan and Kumar |^ first obtained the probability distribution of the reflection coef- 
ficient for a long system with randomness and gain, which is given by: 
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where x = and q = ^o/^i- We numerically calculated the P{x) (or P{R)) for both 
models in cases that q changes drastically. Our numerical features of P{x) {P{R)) give some 
interesting new results. According to Eq.(PBD, the maximum probability of P{x) (or P{R)) 
should appear at 0.5 (or i? = g + 1), and the distribution has a long tail of large x (or R). 

For the many-layered model, when the wavelength is near the band center (A = 360nm), 
we have that L^/^i — 7 for different gains from Eq.(|l^). When q at the range of 0.01 to 50, 
and so < L'f/^i, the P{x) ( P{R)) behaves as predicted by Eg . (^) . When q is at the 
range of 50 to 500, the position of the maximum P{x) {P{R)) begins to shift left away from 
the point 0.5 (or g + 1 ). When q increase further, such as to become close to a thousand, 
the maximum of P{R) shift to the value of Rq, which is the saturated value of reflection 
for the periodic system, and P{R) begin to change its shape into a delta function and the 
long tail disappears. This process is clearly shown in Fig. 7. It is reasonable to assume that 
when q is very large, the system is similar to a periodic system , so the P{R) changes to a 
delta function which is the distribution of the reflection coefficient of the periodic system. 

For the tight-binding model, when frequency is at band center {E = 0), L^/^i is not 
a constant as rj changes. It has a range from 5 to 7. When q at the range of 0.01 to 
20(^0 < Lf/^i), the P{x) ( or P{R) ) behaves as predicted by Eq.dl). When q is at the 
range of 20-400, then the position of the maximum P{x) ( P{R) ) begins to shift left away 
from theoretical value 0.5 ( or g + 1). When q is larger than 400, P{R) develops two peaks 
, one peak evolves from the original peak, the other one emerges at Rq. When q is even 
larger, such as thousands, then the original peak goes down and disappears, and the new 
peak become higher at Rq. At the same time the long tail disappears, the P{R) also changes 
to a delta function at the position of Rq. All these changes are shown clearly in Fig. 8. In 
Ref. |]TB[, they also got two peaks for P{R), but they did not explain that the new peak is 
due to the periodic background of the system and that the delta function is at the position 
of Ro, the saturated value of periodic system. When q is very small, we obtain that the 

14 



P{x) ( or -P(-R)) is almost same as the one predicted by Eg. p6|) , quite different from results 
of [jl3| . We think that this difference is due to the fact that they have not renormalized their 
numerical results. 

In summary, from our numerical results, we get the general behavior of P{x) ( or P{R) 
) for both models. When < -^c^/^b the P{x) ( or P{R) ) is same as the theoretically 
predicted one by Eq.(|26|). When is bigger than we must think about the effect of 

the periodic background and if is really very large, the periodic background will dominate 
the behavior of P{R). We also find that at the band edge wavelength for many-layered model 
(A=470nm) or at the band edge energy for the tight-binding model {E = 1.8), the results of 
P{x) ( or P{R) ) are always different from the predictions given by Eq(p6D. We know that 
at the band edge the effects of coherent scattering will be very strong and make the long 
paths of wave propagation more important . So we think it is this coherent scattering effect 
which make the P{x) ( or P{R) ) different from the prediction of Eq(^) which is obtained 
for homogeneous random systems. 

VI. CONCLUSION 

In conclusion, we have studied the transmission and reflection coefficients in periodic or 
periodically correlated random systems with homogeneous gain. Theoretically, for periodic 
systems we predicted the behaviors of transmission and reflection coefficients, such as the 
slopes of long systems and of short systems and critical length by the transmission matrix 
method. For random systems, first the Zhang's formula of the localization length for long 
systems is checked. We find that only at the band edge and with very strong gain and 
strong disorder, there is obvious deviation from the theoretical prediction of the localization 
length with gain. For short systems, our numerical results show that our generalization 
of the formula of absorbing system is correct for amplifying systems. According to this 
generalization we can predict the behaviors of average value of logarithm of transmission 
coefficient < InT > from the value of 1/^', such as: if it is positive then the < InT > will 
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increase from origin at slope 1/^' and generate a peak at Lc and then start to decrease at 
slope —1/^, if it is negative then the < InT > will decrease monotonically and has a turning 
point at Lc with the slope change from -1/|^'| to — 

To explain the new behavior of the critical length Lc which we got from our numerical 
results, we compare the Letokhov theory with the Lamb theory and give a general expression 
for the critical length considering both the effects of localization and periodic background. 
With this comparison, we also construct the relation of the quality factor Q of a random 
system with the localization length ^. 

We also study the probability distribution of the reflection coefficient P{R) of random 
systems with gain. We find some new behaviors of P{R) and give the criteria for the range of 
validity of the different behaviors and explain it by the influence of the periodic background 
too. The study of wave propagation in amplifying random system is a new and challenging 
topic. There are still a lot things to be done. 
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FIGURES 

FIG. 1. The logarithm of the transmission coefficient T versus (L — L°)/^i, where L° is the crit- 
ical length and is the gain length of periodic systems, (a) For the periodic many-layered model, 
(i), (n) and {in) are the values at three representable wavelength A= 360 nm(band center), 420 
nm(general) and 470 nm(band edge) respectively. The different symbols represent values obtained 
from different gains, e" = -0.001, -0.002, -0.005, -0.001, -0.1. (b) For the periodic tight-binding 
model with E=0 for different gains r?=0.01, 0.05, 0.1, 0.2, 0.5. 

FIG. 2. The logarithm of the reflection coefficient R versus L for the periodic many-layered 
model, (a), (b) and (c) are values of three representative wavelengths A= 360 nm, 420 nm and 470 
nm respectively. From right to left, the numbers on the peaks are the values of Lj?, corresponding 
to different gains e" = -0.001, -0.002, -0.005, -0.01, -0.1. Notice that the saturated value of R is 
independent of e" for the three wavelengths studied. 

FIG. 3. The logarithm of the reflection coefficient R versus L for periodic tight-binding model, 
(a), (b) and (c) are values of three representative energy E =0, 1, 1.8. From right to left, the 
numbers on the peaks are the values of corresponding to different gains ry =0.01, 0.05, 0.1, 0.2, 
0.5. Notice that the saturated value of R for each E depends on rj. 

FIG. 4. The average values of logarithm of T versus L/^. The results the random many-layered 
model, with A=360 nm and W=0.2, are shown by solid lines. Lines from lower to higher correspond 
to different gains e" = 0.0005, 0.001, 0.005, 0.01. When e" is equal to 0.001, l/^i is almost same as 
1/^0) so < In T > is almost horizontal for small L, as shown by the wide solid line. For e" > 0.001, 
1/^1 > 1/Co) and for e" < 0.001, l/^i < l/^o- Results for the random tight-binding model, with 
E = and W = 1, are shown by dashed lines. Lines from lower to higher correspond to r/= 
0.01, 0.02, 0.08, 0.3. When 77 is equal to 0.02, l/^i is almost same as l/^o, so < InT > is almost 
horizontal for small L, as shown by the wide dashed line. For r) > 0.02, l/^i > l/^o, and for 
r? < 0.02, 1/a < 1/^0- 



19 



FIG. 5. 1/^ versus to l/^o + l/Ci; where ^ is localization length for a system with gain, is 
the localization length of a system with disorder but with zero gain, and .^i is the gain length, (a) 
For the random many-layered model, empty symbols are of wavelength A=360 nm, filled symbols 
are of A=470 nm. Different symbols represent different sets of parameters of disorder and gain, 
W = 0.05,0.1,0.2,0.5; e"=0.001, 0.002, 0.005, 0.01, 0.05, 0.1. (b) For the random tight-binding 
model, empty symbols are of energy = 0, filled symbols are of energy E = 1.8. Different symbols 
represent different sets of parameters of disorder and gain, VF=0.5, 0.8, 1, 1.5, 2, 3, 5; r/=0.01, 
0.05, 0.1, 0.3. 

FIG. 6. The critical length is plotted versus different random strengths W . (a) For the 
random many-layered model with A=360nm, the dashed line and darkened line are the values 
obtained according to Zhang's formula and Eq.(p^) respectively, (b) For the random tight-binding 
model with E={) and 77=0. 01, the dashed line and darkened line are the values obtained according 
to Zhang's formula and Eq.(|^) respectively. In both cases Q = and a = 0.7. 

FIG. 7. Probability distribution of the reflection coefficient P(x) versus x of the random 
many-layered system with gain at A = 300 nm for q = 1.1 and 17.7 (a); q = 0.163 (b); q = 451 (c), 
where x = {R — l)/g and q = ^0/^1 • The solid curve given by the solid line in (a), (b) and (c) is 
the analytical result of Eq.(26). In (d), P{R) versus R is plotted for two values of q , q = 1800(low 
one) and 7200(high one). Notice that P{R) approaches a delta- function distribution at Rq when 
q = 7200. 

FIG. 8. Probability distribution of the reflection coefficient P(x) versus x of the random 
tight-binding system with gain at = for q = 1.0 (a); q = 132 (b); q = 525 (c), where 
X = {R — l)/q and q = ^o/Ci- The solid curve given by the solid line in (a), (b) and (c) is the 
analytical result of Eq.(26). In (d), P{R) versus R is plotted for q = 5.25 x 10^. P{R) approaches 
a delta-function distribution at Rq when q = 5.25 x 10^. 
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